Analytical calculation of fragmentation transitions in adaptive networks 
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In adaptive networks fragmentation transitions have been observed in which the network breaks 
into disconnected components. We present an analytical approach for calculating the transition 
point in general adaptive network models. Using the example of an adaptive voter model, we 
demonstrate that the proposed approach yields good agreement with numerical results. 
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In the past decade networks have proven to offer a 
metaphor for describing and analyzing complex systems 
in different fields with applications ranging from techno- 
logical to biological and social systems [i|-lJl. By con- 
ceptually reducing a complex system to a set of discrete 
nodes connected by links a simplification is achieved that 
often enables deep insights in the structure and dynamics 
of the system. 

In network physics, dynamics can refer to two different 
concepts. First, the dynamics OF networks describes the 
temporal evolution of the topology, the specific pattern 
of nodes and links. Second, the dynamics ON networks 
refers to the evolution of internal properties in the net- 
work nodes, which are coupled according to the (typi- 
cally static) topology. Systems combining both types of 
dynamics are called adaptive or coevolutionary networks 

Si- 

Adaptive networks are presently studied in several dif- 
ferent disciplines because they appear in many differ- 
ent applications Q and exhibit some unique phenomena, 
including robust self-organization to critical states Q, 
emergence of distinct classes of nodes from an initially 
homogeneous population [9|, formation of complex hi- 
erarchical topologies [IS El and complex network-level 
dynamics and phase transitions [l2j . 

One exciting recent discovery in the field of adaptive 
networks is the existence of a generic scenario for frag- 
mentation transitions [isj . In many adaptive networks, 
dynamics is contingent on the presence of so-called ac- 
tive links that connect nodes in different states. When 
no such links are left, dynamics stops and an absorbing 
state is reached. Absorbing states are therefore encoun- 
tered either when the network becomes polarized so that 
all nodes are in the same state or when the network frag- 
ments such that nodes separate into disconnected com- 
ponents, which are internally state-uniform. In the fol- 
lowing we denote the phase transition that separates the 
fragmented phase from the polarized phase as the frag- 
mentation transition. 

Fragmentation transitions are frequently observed in 
simulations [l3l-[l7l. However, for a detailed understand- 
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ing of dynamics in more complex future models analyti- 
cal approaches to the fragmentation transition will be in- 
strumental. Existing approaches that faithfully capture 
other transitions |12l . ll7H23l | yield only very rough ap- 
proximations for the fragmentation threshold. A paradig- 
matic example is provided by the adaptive voter model 
that describes o pini on formation in networked popula- 
tions [l3 - [l7l [22l. l23l| . In this model mean-field and pair 
approximations overestimate the fragmentation thresh- 
old by 150-200% [l^,!!!,!!- 

In the present paper we propose a simple analytical ap- 
proach for determining fragmentation thresholds in adap- 
tive networks. This approach takes inspiration from per- 
colation arguments that are commonly used in the com- 
putation of the epidemic threshold of infectious diseases. 
Below we start by providing a brief overview of the adap- 
tive voter model, which is subsequently used as an il- 
lustrative example. We then propose a first approach 
for computing the fragmentation threshold in networks. 
This approach relies on the assumption that the degree 
distribution, i.e. the distribution of the number of links 
connecting to the network nodes, is sufficiently narrow. 
Finally, we extend this approach to arbitrary degree dis- 
tributions. 

The voter model was originally proposed as a concep- 
tual model of opinion formation [2J] . The model consid- 
ers a network in which the nodes represent social agents 
and the links represent social contacts. Each agent can 
hold either of two symmetric opinions. The opinions 
are updated continuously by either a) selecting a ran- 
dom agent and letting it adopt the opinion of a randomly 
chosen neighbor (direct voter model), b) selecting a ran- 
dom agent and copying its opinion to a randomly chosen 
neighbor (reverse voter model) or c) selecting a random 
link and letting one of the linked agents adopt the others 
opinion (link-update voter model). In the non-adaptive 
variants of the voter model the network topology remains 
static. If the network is connected the dynamics therefore 
continues until eventually an absorbing consensus state 
is reached in which all nodes hold the same opinion. 

Adaptive extensions of the voter model take an addi- 
tional process into account. At a certain rate agents that 
are connected to agents of different opinion break the re- 
spective link and connect to a randomly chosen agent of 
their own opinion. This rewiring can lead to a fragmen- 
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FIG. 1: (Color online) Illustration of the evolution of active 
links in a degree regular network with degree fe = 3. Agents 
are depicted as nodes which are open or solid depending on 
their opinion. Shown is the network in the neighborhood of 
an active link connecting nodes of different opinions. Arrows 
correpond to dynamical updates and are labeled with the cor- 
responding transition rate. Depending on the parameters the 
updates lead to proliferation or decline of isolated active links 
(encircled dotted) and 2- fan motifs (encircled dashed). 



tation of the network, such that both opinions survive in 
disconnected network components that are internally in 
consensus [13: J6, 25, 26]. Below, we consider an adaptive 
voter model with link update in which rewiring events oc- 
cur at a rate p and opinion updates occur at a rate 1 — p. 
The challenge in this model is to compute the fragmen- 
tation threshold, i.e. the critical value of p at which the 
active links disappear. 

To simplify the discussion below we define the notion of 
a g-fan motif: a subgraph consisting of one node holding 
a given opinion and q neighbors of the node holding the 
other opinion. The node holding the solitary opinion is 
denoted as the base node of the fan, whereas the other 
nodes are denoted as the fringe nodes. For instance a 
4-fan contains a base node and 4 fringe nodes, which are 
connected to the base node by active links. 

For computing the fragmentation transition we con- 
sider a single active link connecting two almost discon- 
nected clusters of different opinions. We then study the 
dynamics in the vincity of this link and compute the net 
balance of active links. If this balance is negative then 
the number of active links declines exponentially leading 
eventually to the fragmented absorbing state. If the bal- 
ance is positive then active links proliferate preventing 
the network from reaching the fragmented state. 

For simplicity let us first assume that the network is 
degree regular, so that every node has excactly the same 
number of neighbors k. For the special case of fc = 3 
the dynamics of active links is illustrated in Fig. [TJ We 
start from a single active link (left half of figure). In 
the next update that affects this link, the link will be 
rewired with probability p becoming inert (not shown), 
with the complementary probability 1 —p one of the nodes 
connected by this link adopts the other's opinion. In the 
adoption event the original active link becomes inert, but 
the k~l other connections of the adopting agent become 
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FIG. 2: Illustration of the transitions for the g-fan motif {q}, 
in degree regular networks with degree k. Motifs and transi- 
tions are depicted as in Fig. 1. Additionally, the correspond- 
ing systems of differential equations and the corresponding 
Jacobians are shown (see text). 



active. 

The active links that are created by an adoption event 
form a fc — 1-fan. Because these links can become inert 
simultaneously if the node at the base of the fan changes 
its opinion, the links of the fan cannot be considered 
independent. We therefore continue by studying how the 
next update alters the fc — 1-fan motif (Fig. [T] right half). 
If the update is a rewiring event (probability p) then it 
decreases the width of the fan turning the fc — 1-fan into 
the fc — 2-fan. If the update is an opinion adoption event 
(probability \ — p) then there are two possible scenarios 
occuring with equal probability. In the first scenario the 
node at the base of the fan changes its opinion. In this 
case the whole fc — 1-fan becomes inert, but one new 
active link is formed at the base of the fan. In the second 
scenario one of the fringe nodes of the fan adopts the 
base node's opinion, in this case the width of the fan is 
reduced by one, but an additional fc — 1-fan is activated. 

For fc = 3, a fc — 2-fan is just a single active link. We 
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thus obtain a closed system of transitions. Summarizing 
the results of Fig. [U we find that an update affecting a 
single active link, always deactivates the link and acti- 
vates a 2-fan with probability 1 — p. An update affect- 
ing a 2-fan always turns the 2-fan into a single active 
link and additionally creates one 2-fan with probability 
(1 —p)/2 (Fig. m top panel). These observations can be 
summarized further in a system of differential equations 
decribing the motifs 



{1} 



-{l} + 2{2} 



{2} = -2{2} + (l-p){l} + 2.i(l-p){2}, 



(1) 



where {q} denotes the density of q-fans. Because we as- 
sumed that the density of active motifs is low the equa- 
tions of motion are linear and can therefore be solved 
straightforwardly. In particular the fragmented state 
{1} = {2} = is stable and globally attractive if all 
eigenvalues of the corresponding Jacobian 
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have negative real part. The phase transition of the net- 
work thus corresponds to a bifurcation of the low dimen- 
sional model in which an eigenvalue acquires a positive 
real part. For fc = 3 this transition occurs at p = 1/3. 

Degree regular networks with fc > 3 can be treated 
analogously to the fc = 3 example. The corresponding 
equations for the k ~ A and the generalization to arbi- 
trary fc is shown in Fig. [2l An update affecting an active 
link deactivates the link and activates a fc — 1-fan with 
probability 1 — p. An update affecting a g-fan either a) 
deactivates the fan and activates a single link (probabil- 
ity (1 — p)/2); b) decreases the width of the fan by 1, 
turning the g-fan into a, q — 1-fan (probability p); or c) 
decreases the width of the fan by one and activates a new 
fc — 1-fan (probability (1— p)/2). This leads to asystemof 
fc — 1 linear differential equations and the corresponding 
Jacobian. From this Jacobian the transition points can 
then be computed either by computation of eigenvalues 
or by direct analytical construction of testfunctions as 
described in . 

The assumption of degree regularity is questionable be- 
cause even if present in the initial network, degree reg- 
ularity is destroyed by adaptive rewiring. Nevertheless, 
the comparison in Fig. [3] shows that the procedure pro- 
posed above yields significantly better results than a com- 
mon pair approximation. For obtaining a closer match 
with simulation results we propose a refined procedure 
that takes the effect of broader degree distributions into 
account. Instead of the fan motif used above, we now 
consider spider motifs, containing one central base node 
which is connected to m nodes of it's own opinion and 
I nodes of opposing opinion. The {m, Z}-spidcr thus has 
m inert links and I active links, leading to a total degree 
d ^ m + 1. As before we only consider active spiders, i.e. 
I > 0. 

The effects of updates on a spider motif are shown 
schematically in Fig. [D In a rewiring event either the 
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FIG. 3: (Color online) Fragmentation threshold in the adap- 
tive voter model. Shown are numerical results from agent 
based simulation (black dots), pair-approximation (dashed 
line), the analytical approach proposed here (crosses), and 
its refined variant (circles). The proposed approaches yield a 
better match with the numerical results than the established 
procedure. Numerical simulations used = 10^ nodes. In 
the refined approximation a cut-ofi^ of fcmax = 30 was used. 



rewired link is kept by the fringe node and the {m, /}- 
spider is turned into a {m, Z — l}-spider or the rewired 
link is kept by the base node turning the {m, Z}-spider 
into a {m + 1,1 — l}-spider. In an opinion adoption event 
either the base node is convinced, which turns all ac- 
tive links into inert links and vice versa, leading to a 
{Z, m}-spider, or one of the fringe nodes is convinced by 
the base node giving rise to a new {1, g}-spider while in 
the remaining spider one active link turns into an inert 
link. In the latter case g, the number of active links of 
the newly activated spider is given by the next-neighbor 
degree distribution P{g + 1). In the adaptive network 
this distribution is reshaped by the rewiring events. The 
exact distribution that is observed in the long-term be- 
havior is therefore not known. However, the relatively 
good results obtained already with the degree-regular 
approximation suggest that a precise knowledge of the 
distribution is not necessary. In the present example we 
therefore use a Poisson distribution with mean degree fc, 
P{x) = e~'^k^/xl, which adequately captures the effect 
of random rewiring. The transition rules can then be 
summarized as 
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FIG. 4: Example for the transition probabilities of a {3, 4}- 
spider motif. Diagrams and arrows represent motifs and tran- 
sitions in analogy to Fig. 2. The dashed links indicate a num- 
ber of second neighbors which is drawn from a distribution 
(here, Poisson distribution). 



For obtaining a Jacobian of finite size we consider only 
motifs with k ^ m + I < fcmax- Following the same steps 
as above this leads to a (D x D) Jacobi matrix, where 
D = /cmax(fcmax + 1)/2. Because of its complexity we con- 
struct this matrix by computer algebra and compute the 
fragmentation threshold by numerical eigenvalue compu- 
tation. 

The comparison in Fig. [3] shows that the results from 
the refined procedure are in excellent agreement with 
simulation results. In particular the approximation 
works well for low k where pair and mean-field expan- 
sions yield poor results. For high k there is a small dis- 
crepancy, which may be due to long-range correlations or 



the precise shape of the degree distribution. 

In summary we proposed two approaches for comput- 
ing fragmentation thresholds in adaptive networks. The 
simpler approach allows for a quick analytical estimation 
of the threshold with higher accuracy than previously 
proposed approaches. The refined approach yields nu- 
merical predictions with high accuracy. In contrast to 
other techniques the proposed approaches lead to lin- 
ear ODEs because they assume a low density of active 
links. While the linear ODEs can be solved exactly for 
any density of active links, the approximation must fail 
when the system is far from the fragmentation point. 
A difHculty not addressed in the present paper is ob- 
taining the degree distribution for an adaptive network 
without explicit simulation. In principle the relevant in- 
formation for applying the refined procedure proposed 
here could be obtained by third-order moment expan- 
sions |17| . which cannot predict the fragmentation tran- 



sition precisely, but allow relatively faithful estimation of 
the width of the degree distribution. In practice it may 
be simpler to use statistical fits of distributions observed 
in a small set of exploratory simulation runs. Further, 
using a Poisson distribution should yield good results for 
all networks with exponentially decaying degree distribu- 
tions. We therefore believe that both the simple and the 
refined procedure proposed here can be applied with rel- 
ative ease in practice and will be instrumental to explor- 
ing fragmentation transitions in future adaptive network 
models. 

The authors thank G. Demirel and F. Vazquez for in- 
sightful discussions. 
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